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Abstract. A spectral mapping theorem is proved that resolves a key problem in applying 
invariant manifold theorems to nonlinear Schrodinger type equations. The theorem is ap- 
plied to the operator that arises as the linearization of the equation around a standing wave 
solution. We cast the problem in the context of space-dependent nonlinearitics that arise in 
optical waveguide problems. The result is, however, more generally applicable including to 
equations in higher dimensions and even systems. The consequence is that stable, unstable, 
and center manifolds exist in the neighborhood of a (stable or unstable) standing wave, such 
as a waveguide mode, under simple and commonly verifiable spectral conditions. 



1. Main Results 

The local behavior near some distinguished solution, such as a steady state, of an evolution 
equation, can be determined through a decomposition into invariant manifolds, that is, 
stable, unstable and center manifolds. These (locally invariant) manifolds are characterized 
by decay estimates. While the flows on the stable and unstable manifolds are determined by 
exponential decay in forward and backward time respectively, that on the center manifold is 
ambiguous. Nevertheless, a determination of the flow on the center manifold can lead to a 
complete characterization of the local flow and thus this decomposition, when possible, leads 
to a reduction of this problem to one of identifying the flow on the center manifold. 

This strategy has a long history for studying the local behavior near a critical point of an 
ordinary differential equation, or a fixed point of a map, and it has gained momentum in the 
last few decades in the context of nonlinear wave solutions of evolutionary partial differential 
equations. Extending the ideas to partial differential equations has, however, introduced a 
number of new issues. In infinite dimensions, the relation between the linearization and the 
full nonlinear equations is more delicate. This issue, however, turns out to be not so difficult 
for the invariant manifold decomposition and has largely been resolved, see, for instance, [2], 
[3]. A more subtle issue arises at the linear level. All of the known proofs for the existence 
of invariant manifolds are based upon the use of the group (or semigroup) generated by 
the linearization. The hypotheses of the relevant theorems are then formulated in terms of 
estimates on the appropriate projections of these groups onto stable, unstable and center 
subspaces. These amount to spectral estimates that come directly from a determination of 
the spectrum of the group. However, in any actual problem, the information available will, 
at best, be of the spectrum of the infinitesimal generator, that is, the linearized equation 
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and not its solution operator. Relating the spectrum of the infinitesimal generator to that 
of the group is a spectral mapping problem that is often non-trivial. 

In this paper, we resolve this issue for nonlinear Schrodinger equations. We formulate the 
results for the case of space-dependent nonlinearities in arbitrary dimensions. This class of 
equations is motivated by the one space dimension case that appears in the study of optical 
waveguides, see [17], and has attracted the attention of many authors. In particular, there 
is extensive literature on the existence and instability of standing waves, see, for instance, 
[9, 10, 11, 14, 15] and the references therein. In many instances the questions of the existence 
of standing waves and the structure of the spectrum of the linearization of the nonlinear 
equation around the standing wave are well-understood, see [17]. 

In the case considered in this paper, the interesting examples are known to have the 
spectrum of their linearization A enjoying a disjoint decomposition: the essential spectrum is 
positioned on the imaginary axis, and there are several isolated eigenvalues off the imaginary 
axis, see [9, 14, 15]. However, as mentioned above, this spectral information about the 
linearization A is not sufficient to guarantee the existence of invariant manifolds. The general 
theory gives the existence of these manifolds for a semilinear equation with linear part A only 
when the spectrum of the operator e*- 4 , t > 0, rather than that of A, admits a decomposition 
into disjoint components. It is not a priori clear that the spectrum of the operator e*" 4 
is obtained from the spectrum of A by exponentiation. Indeed, in the present case, the 
operator A does not generate a semigroup for which this property is known (such as for 
analytic semigroups). Thus, we prove such a spectral mapping theorem (cf. Theorem 1) in 
this paper. This spectral mapping theorem is derived from a known abstract result in the 
theory of strongly continuous semigroups of linear operators (see Theorem 3). To apply this 
abstract result one needs to prove that the norm of the resolvent of A is bounded along 
vertical lines in the complex plane. The corresponding proof is based on Lemmas 5 and 6. 
The main technical step in the proof of these lemmas concerns a result about the high-energy 
decay of the norm of a Birman-Schwinger-type operator (cf. Proposition 8), a well-known 
device borrowed from quantum mechanics. 

We consider the following Schrodinger equation with space-dependent nonlinearity, 

(1) iu t = Au + f(x, \u\ 2 )u + (3u, u = u(x, t) e C, xeR n , t>0, (3 e R, 
n 

where A = denotes the Laplacian, n G N, and / is real- valued. Rewriting (1) in 

j=i Xj 

terms of its real and imaginary parts, u = v + iw, one obtains 

v t = Aw + f(x, v 2 + w 2 )w + /3w, 
w t = —Av — f(x, v 2 + w 2 )v — (3v. 

A standing wave of frequency f3 for (1) is a time-independent real- valued solution u = u(x) 
of (2). Suppose the standing wave u is given a priori. Consider the linearization of (2) 
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around u (recalling w — 0), 

p t = Aq + f(x,u 2 )q + (3q, 

q t = -Ap - f(x, u 2 )p - 2d 2 f(x, u 2 )u 2 p - ftp, 

where d 2 f(x, y) = f y (x,y). Thus, the linearized stability of the standing wave is determined 
by the operator 



(3) A 



-L R 
Li 



whereL fl = -A-/3 + Qi, L 7 = -A - (3 + Q 2 , 



and the potentials Q\ and Q 2 are explicitly given by the formulas 

Qi(x) = -f(x, u 2 (x)), Q 2 (x) = -f(x, u 2 (x)) - 2d 2 f(x, u 2 (x))u 2 (x). 

We impose the following conditions on f,f3 and the standing wave u (see [9, 14, 15]): 

(HI) / : M. n+1 — > K. is C 3 and all derivatives are bounded on a set of the form K. x U, where 

U is a neighborhood of G MP; 
(H2) f(x, 0) — > exponentially as |x| — > 00; 
(H3) (3 < 0; 

(H4) |"u(x)| — > exponentially as \x\ — > 00. 

As a result, the potentials Qi and Q 2 exponentially decay at infinity. The operator A is 
considered on L 2 (M. n ) © L 2 (R™); the domain T>{— A) is chosen to be the standard Sobolev 

A 



space # 2 (M n ), and the domain of A is then H 2 (R n ) © iJ 2 (IR n ). Note, that 



-A 



generates a strongly continuous group on L 2 (R n ) ©L 2 (M n ). Thus, its bounded perturbation 
A generates a strongly continuous group {e t,A }teM. as well. 

It was proved in [9, Thm. 3.1] that cr css (^4) = {i£ : £ G K, |£| > — In addition, it 
was proved in [9, 14, 15] that, under the above hypotheses, a(A)\a ess (A) consists of finitely 
many eigenvalues, symmetric with respect to both coordinate axes. 

We prove the following result that relates the spectrum of the semigroup {e tA } t > and the 
spectrum of its generator. 

Theorem 1 (Spectral Mapping Theorem). For each n G N one has 

a(e tA )\{0} = for allt>0. 

See [8, 20, 23, 25] for a discussion of the spectral mapping theorems for strongly continuous 
semigroups and examples where the spectral mapping property as in Theorem 1 fails. Since 
the spectral mapping theorem always holds for the point spectrum, Theorem 1 implies, in 
particular, that a css (e t - A ) C T = {\z\ — 1}. A spectral mapping theorem was proved in [18] 
for a nonlinear Schrodinger equation with a specific potential and in the case n — 1. Their 
proof also uses Theorem 3 that is the key to our result. To the best of our knowledge the 
work of Kapitula and Sandstede [18] was the first to use the Gearhart-Greiner-Herbst-Pruss 
Theorem in this context. 
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It follows that there will be only finitely many eigenvalues of e tA off the unit circle and 
therefore general results on the existence of invariant manifolds for semilinear equations can 
be invoked (see, e.g., [3] and compare also with [4] and the literature cited in [19, p. 4]). 

The (local) stable manifold is defined as the set of initial data whose solutions stay in the 
prescribed neighborhood and tend to u exponentially as t — > +oo. The unstable manifold is 
defined analogously but in backward time. The center manifold is complementary to these 
two and contains solutions with neutral decay behavior (although they can decay, they will 
not do so exponentially). In particular, the center manifold contains all solutions that stay 
in the neighborhood in both forward and backward time. For details see, for instance, [3]. 

Concerning the equations under consideration here, we have the following main theorem. 

Theorem 2. Assuming (Hl)-(H^), in a neighborhood of the standing wave solution u of 
(1) there are locally invariant stable, unstable and center manifolds. Moreover the stable 
and unstable manifolds are of (equal) finite dimension and the center manifold is infinite- 
dimensional. 

All of the examples of standing waves considered in [16] and [17] satisfy the hypotheses 
given here and thus enjoy a local decomposition of the flow by invariant manifolds. Some 
of these waveguide modes are stable, while others are unstable. In the unstable cases, the 
above results show that the instabilities are controlled by finite-dimensional (mostly, just 
one-dimensional) unstable manifolds. A natural question is whether the waveguide modes 
are stable relative to the flow on the center manifold. Such a result was shown for the case 
of nonlinear Klein-Gordon equations in [4] using an energy argument. Whether such an 
argument will work for nonlinear Schrodinger equations is open. It is more than of academic 
interest, as stability on the center manifold has the consequence that the center manifold is 
unique, see [3], and armed with such a result, a complete description of the local flow can be 
legitimately claimed. Cases of standing waves in higher dimensions are given in [15]. Some 
of these are unstable and the above considerations again apply. 

We also wish to stress that the spectral mapping theorem developed here is not restricted 
to a single equation. Indeed, the results formulated here are easily adaptable to systems 
of nonlinear Schrodinger equations. This is particularly important as such systems arise 
in, among other problems, second harmonic generation in waveguides and wave-division 
multiplexing in optical fibers. The case of systems is considered in Section 4. 

In the next section, we give the basic set-up that will be used and formulate the necessary 
lemmas for proving Theorem 1. The proofs are given in Section 3. 

2. Basic lemmas 

Fix an a G K\{0} such that the line {£ = a + it : r G M} does not intersect a (A). To 
prove Theorem 1, one needs to show that a(e tA ) does not intersect the circle with radius 
e ta centered at the origin. We will use the following abstract result known as the Gearhart- 
Greiner-Herbst-Priiss theorem. 
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Theorem 3. (see, e.g., [20, p. 95].) Let A be a generator of a strongly continuous semigroup 
on a complex Hilbert space. Assume 

a(A) n {£ = a + ir : r G K} = 0, a G IR\{0}. 

Then a{e tA ) fl {|z| = e* a } = /or a// 1 > 0, if and only if the norm of the resolvent of A is 
bounded along the line {£ = a + ir : r G IR} , £/ia£ is, 

(4) sup || (a + ir - A)' 1 ]] < oo. 

Therefore, to prove Theorem 1, one needs to show that || (a + ir — »4) _1 || remains bounded 
as |r| — > oo for the operator A defined in (3). 

We denote D = —A — f3 and recall that (3 < by (H3). Moreover, we have 

<t(D) = <t(-A)-P = [-0,oo). 

We note, that D 2 with domain H^iW 1 ) is a self-adjoint operator. Thus, for £ = a + ir with 
r^O one has — £ 2 ^ cr(D 2 ). Moreover, we write 

-L T t { -D £ I + I -g 2 o 



(5) 



£ D \ 7 Q 1 
-D i ) \-Q 2 



where, by a direct computation with operator-valued matrices, 

i d \ 1 _ f ^e+D 2 }- 1 -ie+D 2 ]-^ 

Lemma 4. Fort; = a + ir, a G M\{0}, reK, the norm of the operator (6) remains bounded 



(6) 



as r — > oo. 



The elementary proof of this lemma is given in the next section. 
Next, we denote 

m T(n _f t DyW Q 1 \_\[e + D 2 ]- l DQ 2 ^ + £> 2 ]- 1 Q 1 

(7j V —D Cj [-Q2 J " [ -tit 2 + DT'Qz [£ 2 + D 2 ]~ 1 DQ 1 

The main step in the proof of Theorem 1 is contained in the next two lemmas. They imply 
that the norm of the operator (i+T(£)) _1 is bounded as |r| — > oo. Assume Q is a real-valued 
potential exponentially decaying at infinity and let £ = a + ir. 

Lemma 5. (a) For n = l, ||Q[£ 2 + -D 2 ] _1 D|| -> as |r| -> oo; 
(b) For n > 1, ||£[£ 2 + D 2 ]^ 1 ^!! -> as |r| -> oo. 

Lemma 6. For eac/i n > 2 there exists a r > snc/i t/iat /or |r| > r £ae operator I + [£ 2 + 
F> 2 ] -1 F>(5 aas a bounded inverse. Moreover, 

(8) sup HfJ+^ + F* 2 ]- 1 ^]- 1 !! <oo. 

|t|>t 
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The proofs of Lemmas 5 and 6 are given in the next section. We proceed finishing the 
proof of Theorem 1. 

Proof of Theorem 1. In the case n — 1, by passing to the adjoint operator of £[£ 2 + D 2 ]~ 1 Qi 
and [£ 2 + D 2 ]~ 1 DQ 2 , parts (a) and (b) of Lemma 5 imply that the norm of each of the four 
block-operators in the right-hand side of (7) is strictly less than 1 for |r| sufficiently large. 
Thus, ||T(0|| < 1. By (5), one infers 



(a + ir-AY 



-D 



D 



< 



i-imoii 



-D £ 



Using Lemma 4, one obtains (4) and hence Theorem 3 implies the result. 
For n > 2, we write I + T(f ) = I + M(£) + JV(f ), where 



M(0 
N(0 



\e + D 2 Y 1 DQ 2 

[e + D^DQ, 

o tte+^Qi 

i\e + D^Q 2 



as 



oo. 



Lemma 6 yields that 
-l i 



By part (b) of Lemma 5, one has ||iV(£)|| 

I + M(£) is invertible for \t\ > r , with K := sup| r |> ro ||(J + M(^)) _1 || < oo. Choosing 
r o — r o such that lf||jV(£)|| < 1 for |r| > Tq, the norm of the operator (7 + T(^)) -1 = 
[/+(/ + M(0) _1 ^V(0] _1 (^ + M^))- 1 is bounded for |r| > and, as before, Theorem 3 
implies the result. □ 



3. Proofs 

In this section we give the proofs of Lemmas 4-6. The proof of Lemma 5 is based on a 
direct estimate of the trace norms. We will give two proofs of Lemma 6. The first proof is 
applicable to all n > 2 and uses estimates for the norm of the resolvent on weighted spaces 
of L 2 -functions. The second proof works for the cases n = 2 and n = 3 and is based on 
explicit estimates for the integral kernel of the resolvent of the Laplacian. 

The main tool in the proof of parts (a) and (b) of Lemma 5 is the following well-known 
result. Denote by J q the set of operators A e £(L 2 (M n )), n > 1 such that \\A\\ Jq = 
(ti(\A\ q )) 1 / q < oo, q > 1, where tr(-) denotes the trace of operators in L 2 (M n ). We recall that 
\\A\\ < \\A\\ Jq for all q>\. 

Theorem 7. (see, e.g., [21, Theorem XI.20].) Suppose 2 < q < oo and let f, g e L q (W n ). 
Then f(x)g(-iW) G J q and 

||/(^HV)b 9 <(27r)-"/ 2 ||/|| L9 ||^|| L ,. 
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We will apply Theorem 7 to the exponentially decaying Q = f E L q , q > 1 and an 
appropriate choice of g. 

Proof of Lemma 5. First, let g(x) = (\x\ 2 — /3)(£ 2 + (\x\ 2 — /3) 2 ) -1 . Then, for n > 1, one has 
g(-iW) = [£ 2 + (-A - l3) 2 ]-\-A -13)= g{-id/dx) = [£ 2 + D 2 ]~ L D. For r = \x\ one infers, 

POO 

/ r n-l( r 2 _ /^[(^ _ r 2 + (r 2 _ ^2)2 + 4^2]-^ 
JO 

(r 2 -(3) 2 -(r 2 -a 2 )^ 2 



Li 



(2\ar)r f 
Jo 



r n - l (r 2 -py 



2\ar\ 



The change of variable y = ((r 2 — (3) 2 — (r 2 — a 2 ))/(2\ar\) shows that 



h = k(r) 



{y 2 + i)- q,2 [-,y + i 



e" 1 

2\ar\ 



1 



2-1 

2 



where 



and 



, n 1 „, O On g-1 = 



T 



7 y + i + 



+ 1 







-9/2 



dr. 



(r 2 — a 2 ) 2 



n-2 
2 



fc(r) = ^(2|ar|) 1 -"(r 2 - a 2 )^+^ = 0(\r\^) as |r| -> oo 



r' = (r 2 - a 2 )/(2|ar|) = 0(|r|) as |r| -> oo. 
We recall, that (3 < by (H3). Thus, 

(9) \\g\\l q <c\r\^I(-r',oo), 
denoting 



I(a, b)= / (y 2 + l)^ 2 l-y + 1 



q—l I re — 2 
2 4 



Next, one observes that Irl^/fY', oo) — > as |t| — > oo for g > max{l,|}. Indeed, for 
y > t' one has 

/CO 2g + n-4 

( y 2)-*/2 ( 2 £) 4 ^ = c | T | 2 ^r' 1 -^0as|r|^oo. 

n — 2 

To prove part (a) of Lemma 5 for n = 1 we remark that \t\^I(—t', t') — > as |r| — > oo. 
We stress, that this assertion does not hold for n > 2; that is why one needs Lemma 6. 
Indeed, assuming n = 1, for y < r' one concludes 

/•oo 

|r|-5/(-r', r) < 2 ^- 3 )/ 4 |r|-5 / + iy*/* dy _> as |r| -> oo. 



Therefore, for n = 1 one infers H^Hl? as |r| — > oo and, using Theorem 7, 

ligr + D 2 ]- 1 ^!! < \\Q(x)g(-iV)\\ Jq < c\\Q\\ Lq \\g\\ Lq ^0 as |r| - oo. 
Thus, part (a) of Lemma 5 is proved. 
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To prove part (b) of Lemma 5, let g(x) = £(£ + (\x\ 2 — j3) ) . Then, for n > 1, one 



obtains g(— iV) = £[£ 2 + (—A — f3) 2 } l . For r = \x\ one concludes as above, 







V 2 — g 2 ); 



n-2 
2 



where 



Z(r) = i(2|ar|) 1 ^(r 2 - a 2 ) 2 ^ = 0{\t\^) as |r| - oo. 



Since (3 < one has 
denoting 



\g\\ Lq < c\t\ 2 J(-t ,oo), 



J(a, 6) 



-9/2 



7 1/ + 1 



Next, fix g > max{l, (n + 2)/2}. We claim that |r| " ^ J(r', oo) — > as |r| — > oo. Indeed, 
for y > r' one concludes, 



n — 4 
4 



oh/. 



n — 2q . 

T\ 2 J(t , oo) < \T 



n — 2q 



(y 



2,-q/2 ( V_ 



4 



dy = c\t\ 2 



as |t| — > oo. 



On the other hand, for y < r' one has 



i n — 2q , , n — 4 . . n — '2q 

\t\ — J(-t',t') < 2~ \t\ — 



/oo 
(y 2 + ir 9/2 cfo/ -> as |t| -> oo. 



□ 



Therefore, \\g\\Li — ► as |r| — > oo and, using Theorem 7, one obtains for n > 1, 

||Q£[£ 2 + (-A - /3) 2 ]- 1 !! < ||g( x )^(-iV)|U 9 < cIIQIIwII^Hl, - as |r| - oo. 

Thus, part (b) of Lemma 5 is proved. 

The proof of Lemma 6 is based on the following proposition. Denote 

(10) \Q\^ 2 (x) = \Q(x)\^, Q^(x) = \Q\^(x) sgn[Q(x)], 
so that Q = Q l / 2 \Q l ' 2 \. 

Proposition 8. Assume Q is a real-valued potential exponentially decaying at infinity. Then 
for n > 2 and Im(cj) > one infers 

(11) IHQr^-A-^-^H^^^O as |Re(cu 2 )| ^oo. 

We postpone the proof of the proposition and proceed with the proof of Lemma 6. 

Proof of Lemma 6. Recall the following elementary fact: If A and B are bounded operators, 
then / + AB is invertible provided / + BA is invertible; moreover, 

(I + AB)- 1 = I-A(I + BA)- 1 B. 
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Thus, using Q = Q^IQI 1 / 2 and letting A = [£ 2 + D^DQ 1 ' 2 and B = \Q\ 1/2 , the desired 
inequality (8) is implied by the following claim 

(12) IIIQI^ + D 2 ]- 1 ^ 1 ^! ^Oas |r|->oo. 

Indeed, if (12) holds, then, for some To > 0, one has \\BA\\ < 1 and sup| r | >ro || (/ + .£? A) < 
oo. To prove sup| T | >To \\A\\ < oo we use the identity 

(13) [e + D 2 ]~ 1 D = \[D + ifl" 1 + \[D - it]' 1 - 
Since D is self-adjoint, 

|| [D ± = 1/| Im(±*OI = 1/M < oo. 

This and \\B\\ = |||Q| 1/2 ||oo < oo implies sup| T | >ro || (/ + AB) ^\\ < oo, which is the desired 
relation (8). 

To prove the claim (12) one recalls that D = —A — (3 and similarly to (13) one obtains 

(14) \Q\ l / 2 \e + D 2 ]~ l DQ l l 2 = i|Q| 1/2 [-A - ^Q 1 ' 2 + i|Q| 1/2 [-A - ^Q 1 ' 2 . 

Here ui\ and u 2 are chosen such that Im(a;i) > 0, Im(aj 2 ) > and uj\ = f3 — i£, lo\ = f3 + i£. 
Recall that £ = a + ir. Thus, if |r| — > oo then | Re(co> 2 )| — > oo for j = 1, 2. An application of 
Proposition 8 to each summand in the right-hand side of (14) then proves the claim (12). □ 

We will give two proofs of Proposition 8. The first proof is based on an estimate for the 
norm of the resolvent for the Laplacian acting between weighted L 2 -spaces. This estimate 
is given in Lemma 9 below. In fact, Lemma 9 is just a minor refinement of Lemma XIII. 8. 5 
in [22]. Results of this type go back to Agmon [1], Ikebe and Saito [12], and others (cf. the 
discussion in [22, p. 345-347]). Much more sophisticated results of this type can be found in 
the work by Jensen [13] and the bibliography therein. 

For n > 2 and s > 1/2 let p s {x) = (1 + \x\ 2 ) s / 2 , for x G R n , and consider the weighted 
L 2 -spaces 

L 2 s (R n ) = {f--\\f\\s--= ||p s /IU W <oo} 

and 

L 2 S (M") = {/: ||/||_,:= || P ;V||l W < oo}. 

Note that L 2 {R n ) ^ L 2 {R n ) ^ L 2 _ s (R n )- Also, let S(R n ) denote the Schwartz class of 
rapidly decaying functions. 

Lemma 9. There exists a constant d = d(n, s) depending only on n and s, such that for all 
A G C with Im(A) ^ and | Re(A)| > 1, and all <p G 5(lR n ) the following estimate holds 

(15) ||^||^<rf|Re(A)r 1 / 2 ||(-A-A) V 9|| s . 
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Proof. For completeness, we briefly sketch a modification of the proof of Lemma XIII. 8. 5 
in [22] to prove inequality (15). An elementary calculation shows that for each A G C with 
Im(A) 7^ the expression 

T 2 ■= mf [\y 2 - A| 2 + \y\ 2 } = inf [y 2 + (1 - 2 Re(A))y + (Re(A)) 2 + (Im(A)) 2 ] 

is equal to Re(A) + (Im(A)) 2 - 1/4 for Re(A) > 1/2 and to |A| 2 for Re(A) < 1/2. Therefore, 
7 < 2|Re(A)|~ 1 / 2 for |Re(A)| > 1. 

For a positive a to be selected below, let pi, a (x) = (1 + alx] 2 ) 1 ^ 2 . Arguing as in the proof 
of Lemma XIII. 8. 5 in [22] (see the corresponding equations (61)-(62) in [22]), one infers 

n 

where d depends only on n and s. Next, pick a < 1 such that da + na l l 2 s < 1/4. Note that 
a depends only on n and s. Moreover, ^(da + na l l 2 s) < 1/2 uniformly for | Re(A)| > 1, 
since 7 < 2| Re(A)|" 1/2 . Since s > 1/2, one obtains 

1 n 

2 WpI,>\\l* < 7(2s + l)(||pl£(-A - \M L 2 + ^ llp^^lUO- 

j'=i 

This and the inequality < < a~ s ^ 2 p~ s show that 

n 

IMU = IIpJVII^ < l|pr>IU» <7(2s + l)a- s/2 (||(-A-A)v9||_ s + ^||9^||_ s )- 

j'=i 

Now (15) follows from the inequalities || • ||_ s < || • \\ s and H^c^H _ s < C|| (—A — A)<p|| s , where 
C is an absolute constant (see Lemma XIII. 8. 4 in [22]). □ 

If Im(A) 7^ 0, then the resolvent (— A — A) -1 is a bounded operator on L 2 (IR n ). Consider its 
restriction R S _ S (X) := (— A — A) -1 | L2 mn) as an operator from L 2 (IR n ) to L 2 _ s (R n ). Inequality 

(15) shows that R s - S (\) is a bounded operator from L 2 s (R n ) to L 2 s (IR n ) and that 

(16) \\Rs,- s mc(LlLl s) < d I Re(A)|~ 1/2 - as |Re(A)| - 00. 

This relation with A = uj 2 will be used in the first proof of Proposition 8. 
First proof of Proposition 8. Define the multiplication operators (M s f)(x) = p s (x)f(x) and 
(M- S f)(x) = p~ 1 (x)f(x). If Q is a real-valued potential exponentially decaying at infinity, 
then using notation (10), one observes that the operators M S Q 1 / 2 and IQ^^Mzl = \Q\ l / 2 M s 
are bounded operators on L 2 (M. n ). On the other hand, M s : £ 2 (M n ) — > L 2 (R n ) and M_ s : 
L 2 _ s (M, n ) — > L 2 (M n ) are isometric isomorphisms. This fact and the identity 

|0| 1/2 (-A-o; 2 )- 1 Q 1/2 /= [IQI 1 ^- 1 ] [M_ s R s ,_ s (u 2 )M^] [M S Q^ 2 ] f 
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for / E L 2 (R n ), yields the following estimate 
IIIQIVVA-u; 2 )- 1 ^ 2 !!^) 

< || |Qr /2 A^r s x || .R s ,_ s (a, 2 ) ll^^g^,) II A^Q 1 ^ | U(i2) . 

Using relation (16) for A = oo 2 one concludes that (11) holds. □ 
The main tool in the second proof of Proposition 8 in the case n = 3 is Theorem 1.23 in 
[24] , inspired by previous results of Zemach and Klein [29] . Assume n = 3 and suppose that 
the potential Q satisfies the following (Rollnik) condition (see [24, p. 3]) 

(17) [ /MfM^oo. 

J J \x-y\ 2 

Note that (17) trivially holds for exponentially decaying continuous Q. Consider on L 2 (R 3 ) 
the operator TZ^ with integral kernel 

e iu\x-y\ 

(18) R„(x,y) = |Q| 1/2 (x) — :Q 1/2 (y), x,yER\x^ y, lm(u) > 0. 

Air\x — y\ 

Theorem 10. ([24, Theorem 1.23].) Assume (17). Then TZ k , k E R, is a Hilbert- Schmidt 
operator in L 2 (M 3 ) with lim|fc|_ >00 ||7?-fe|| = 0, fc 6 M. Moreover, 

(19) tr (n* k n k n* k n k ) -> o as \k\ -> oo, k e r. 

Relation (19) follows from a direct calculation of tr(7l* k 7l k 7l* k 7l k ) , resulting in 




D -ik(\x-y\-\y-z\ + \z-w\-\w-x\) _ |Q(s)||Q(y)||Q(«;)||Q(z)| 



■dxdydzdw, 



(47r) 4 |rc — y\\y — z\\z — w\\w — x 

see [24, p. 24] (and compare also with (22) below), and from the Riemann-Lebesgue lemma. 
The first assertion in Theorem 10 follows from (19) and the estimate ||7?.fc|| 4 = ||7^7?.A;|| 2 < 

\\Kin k \\ 2 j2 = tr(niiz k izin k ). 

In the case n = 2 one replaces the integral kernel e lu) \ x ~y\ / (An\x — y\) in (18) by the Hankel 
function (i/A)H^\uj\x — y\) of order zero and first kind (cf. [6]). We use the following result. 

Theorem 11. ([5], [7, p. 1449].) Assume Q E L 2 (R 2 ) n L 4/3 (M 2 ). Consider on L 2 (R 2 ) the 
operator JC^ with integral kernel 

(20) K u (x,y) = \Q\ l ' 2 {x) l -H^\u\x -y\)Q 1/2 (y), x,yER 2 ,x^ y, lm(u) >0,uj^0. 
Then 

(21) \\!C\\j 2 < c \u;\-^ 2 \\Q\\ L4/3{m , hn(u) > 0, cu + 0. 
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Note that Q G L 2 (R 2 ) n L 4//3 (IR 2 ) for exponentially decaying continuous potentials Q. 

To establish the connection with the discussion in the current paper, we recall that (—A — 
uo 2 )~ l for Im(cj) > is an integral operator with kernel e l ^ x ~ y \ / (4tt\x — y\) in the case n = 3 
and with kernel (i/A)H^\uj\x — y\) in the case n — 2. 

Second proof of Proposition 8. In the case n = 2 we use Theorem 11 for the operator 
K u := \Q\ l / 2 {-A-u 2 )~ l Q 1 / 2 in L 2 (R 2 ) to conclude that |||Q| 1/2 (-A - uo 2 )~ l Q l/2 \\ -> as 
| Re(cj 2 )| — > oo. This proves Proposition 8 for n = 2. 

In the case n = 3 one needs one more calculation. If uj = k + im with m > 0, k e K, then 
the integral kernel S u (x,y) of the Hilbert-Schmidt operator 5 W := |Q|^ 2 (— A — c<j 2 )~ 1 (5 1//2 in 
L 2 (IR 3 ) is given by the formula 

S u {x, y) = R k (x, y)e- m ^-y\ x,yeM. 3 ,x^ y, 

where Rk(x,y) is defined in (18) (with uj = k, Im(c^) = in (18)). 
A direct calculation using m > then shows 

tr(s*s u s*s u ) = [ [ [ [ e -ik(\x-y\-\ y -z\+\z-w\-\w-x\) \Q( x )\\Q(y)\\QM\\Q(z)\ x 

J J J Jri2 \x - y\\y - z\\z - w\\w - x\ 

(22) x {ATxy^e^^y^^-^+^-^dxdydzdw < ti{TV k n k TV k n k ). 

But then relation (19) yields ||<SJ 4 < tr(K* k K k 1l* k 1l k ) -> as -> oo. Similarly, ||<SJ 4 < 
tr( l S* l S wl S* l S UJ ) — > as m — > +oo by the dominated convergence theorem and hence ||»S W || — > 
as |Re(w 2 )| = \k 2 -m 2 \ -> oo. Thus, |||Qj 1 /2[_ A _ ^2]-iqi/2|| ^ as |Re(w 2 )| -> oo, and 
hence Proposition 8 is proved. □ 

Proof of Lemma 4- We recall that £ = a + ir, a e M\{0}, reR. Since D 2 is self-adjoint, 

If I f Q 2 i r 2U/2 

||£(£ 2 + £> 2 ) _1 || < , T 1 ' , = , — r — < c uniformly with respect to \r\. 

11 v | Im(£ 2 )| 2|ar| ~ 1 1 

On the other hand, [£ 2 + D 2 ]~ 1 D = h(D) with ft, defined as h(p) = p(^ 2 + p 2 )^ 1 . Therefore, 

(23) \\(^ + D 2 )- 1 D\\= max \h(p)\ = max |p| [(a 2 - r 2 + p 2 ) 2 + 4a 2 r 2 H. 

pea(D) P>~/3 

To verify that the right-hand side of (23) is bounded as |t| — > oo, one considers two cases. 
First, let p 2 > 2r 2 . Since p 2 - r 2 > \p 2 , 



max \p\ [(a 2 — r 2 + p 2 ) 2 + 4a 2 r 2 ] 2 < max |p| 

|p|>V2|r| ' " 



' 2 , P 2 

a + y 



< V2. 



Second, for p 2 < 2r 2 , one has 



max \p\ [(a 2 - t 2 + p 2 ) 2 + Aa 2 r 2 ] * < max -y^^ < ! 



\p\<V2\t\ " |p|<v^|t| 2|a| |r| y/2\a\ 

This completes the proof of Lemma 4. □ 
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4. Systems of nonlinear Schrodinger equations 



Systems of nonlinear Schrodinger equations arise in many applications of nonlinear optics. 
From the perspective of the current paper, the fundamental issue is the same as for a single 
equation, namely the existence and stability of standing wave solutions. A particularly 
interesting example of multiple pulses was studied recently by Yew [26, 27, 28]. The problem 
of second harmonic generation occurs in a slab waveguide with a quadratically nonlinear 
response. Yew showed that multiple pulses could be generated from the base pulse through 
a resonant homoclinic bifurcation. The details of this process are not pertinent for the current 
work, only the outcome, which is the presence of multiple pulses that Yew has shown to be 
unstable due to real, positive eigenvalues. It can be shown easily that the essential spectrum 
for these pulses remains on the imaginary axis and hence the spectral configuration for the 
linearization at the multiple pulses is analogous to the cases considered above for scalar 
equations. Since the technology of this paper can be applied to this system, Theorem 2 
holds in this case. In the following we show how the above considerations can be adapted to 
this case of systems. 

The equations governing the second harmonic generation problem are a system of coupled 
nonlinear Schrodinger equations of the form 



(24) 



.dw 

i — \- Aw — uw + wv 
at 



io~- 



dv 
di 



+ Av — av H — w 2 



o, 



o, 



where w = w(t,x) and v = v(t,x), x G M n , n > 1, are complex valued functions, and a, 
a, and 9 are positive parameters. For the one-dimensional case n — 1 the questions of the 
existence of standing waves for (24) and, as mentioned above, the structure of the spectrum 
of the linearization around the standing waves are well-understood, see [27, 28]. 

" -L R 
Li 

that has the same structure as in (3), but L R and Li are now 2x2 operator matrices defined 
as follows, 



The linearization at a standing wave u = (</?, iji) is given by the operator A 



(25) 



L R : = 



-A + e- 

-ip/a 



-if 
-A + a) I a 



Li := 



-ip/o 



-if 
-A + a)/cr 



where the functions ip and ip are assumed to be continuous and exponentially decaying at 
infinity. 

Theorem 12. The Spectral Mapping Theorem holds for the group generated on the space 
[L 2 (M n )] 4 by the operator A with L R and Li defined in (25). 

Proof. As above, one needs to show that \\(a + ir — ^4) _1 || remains bounded as 
Denote 

" £>i 



oo. 



D 



Do 



1 a 

where D 1 = —A + 6, D 2 = A + -. 

a a 
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and consider the following matrix potentials exponentially decaying at infinity, 



Qi 



-rj} -if 
-(p/a 



and Q 2 



-<p/a 



With this new notation and ^ = a + ir formula (5) still holds. 

By transposing the second and third rows and columns in the 4x4 matrix 
we remark that this matrix is similar to the block-diagonal matrix with the blocks 

D 2 



—D 



D 



e d 1 

-Dx i 



and 



-D 2 



on the main diagonal and zero remaining entries. Lemma 4, applied to each of the diagonal 
blocks, shows that the norm of the operator (6) in the matrix case remains bounded as 



oo. 



Formula (7) is also valid in the new notation. In particular, one infers, 







By part (b) of Lemma 5, the norm of each entry in the last matrix decays to zero as \r\ — > oo. 
Similarly for £[£ 2 + D 2 ]~ 1 Q 2 , and, as a result, the norms of the off-diagonal entries of T(£), 
defined as in (7), decay to zero as |r| — > oo. 

To handle the diagonal entries of T(£), let Q := [qj,k]\k=i t> e a matrix- valued function 
with exponentially decaying q^. Note that 



[e + D^DQ 



[e+Df]- i D iqil [e+Df]- i D iqi2 

\e + Dl\- x D 2 q 21 [e + Dl]^D 2 q 22 



In the case n — 1, one can apply part (a) of Lemma 5 to each entry of this matrix. This 
shows that ||T(£)|| — > as \t\ — > oo and the proof of Theorem 12 for n — 1 is finished. 

In the case n > 2 we claim that the conclusion of Lemma 6 holds for the matrix- valued case 
(and, hence, Theorem 12 holds). Indeed, to verify (8), one considers the polar decomposition 
Q(x) = \Q(x)\U(x), where \Q(x)\ := [Q*{x)Q{x)\ 1 ' 2 and U(x) is a partial isometry. To 
be consistent with our previous notations in (10), we denote |<5| 1,/2 (^) = ^(a;)! 1 / 2 and 
Q 1 / 2 (x) = \Q(x)\^ 2 U(x), so that Q = \Q\ l / 2 Q 1 / 2 . Similarly to the proof of Lemma 6 in the 
scalar case, we put A = [£ 2 + D 2 ]~ l DQ 1 / 2 and B = {Ql 1 ^ 2 and observe that 

[/ + [£ 2 + D 2 ]~ l DQY 1 = I-A(I + BA)- X B. 

As above, the estimate (8) in the matrix case follows from the relation 

|||g|l/2^2 + D 2]-l £) |g|l/2||_ >0 as | r | _ ^ 

(cf. claim (12)). To verify this relation, we denote |Q| 1//2 (a:) := l a j,k(x)} 2 k=i an d observe 
that the entries of the matrix operator IQ] 1 ^ 2 ^ 2 + D 2 ]^ 1 D\Q\ 1 / 2 are scalar operators of the 
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type 

a hl \e + Df]D iai:k + a h2 \e + Dl]- 1 D 2 a 2tk , j, k = 1,2. 

But the norm of each summand in these entries tends to zero as |r| — > oo. To verify the 
latter fact one can apply identity (14), Proposition 8, and Lemma 9 for the scalar case, and 
the result follows. □ 
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